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In inflationary cosmology, the particles constituting the Universe are created after inflation in 
the process of reheating due to their interaction with the oscillating inflaton field. In the bosonic 
sector, the leading channel of particle production is the non-perturbative regime of parametric 
resonance, preheating, during which bosons are created exponentially fast. Pauli blocking prohibits 
the unbounded creation of fermions. For this reason, it has been silently assumed that the creation 
of fermions can be treated with perturbation theory for the decay of individual inflatons. We 
consider the production of fermions interacting with the coherently oscillating inflatons. We find 
that the actual particle production occurs in a regime of the parametric excitation of fermions, 
leading to preheating of fermions. Fermion preheating differs significantly from the perturbative 
expectation. It turns out that the number density of fermions varies periodically with time. The 
total number of fermions quickly saturates to an average value within a broad range of momenta 
cx q 1//4 , where q is the usual resonance parameter. The resonant excitation of fermions may affect the 
transfer inflaton energy, estimations of the reheating temperature, and the abundance of superheavy 
fermions and gravitinos. Back in the bosonic sector, outside of the parametric resonance bands there 
is an additional effect of parametric excitation of bosons with bounded occupation numberin the 
momentum range oc q 
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In the inflationary scenario, the Universe initially ex- 
pands quasi-exponentially in a vacuum-like state with 
vanishing temperature, entropy, and particle number 
densities. Consider a simple chaotic inflation during 
which all energy is contained in the inflaton field <f> which 
is slowly rolling down to the minimum of its effective po- 
tential V (</>). When inflation ends at <f> ~ M p , the infla- 
ton field begins to oscillate near the minimum of V(<f>) 
with a very large amplitude, <f> ~ JqM p . This scalar field 
can be considered as a coherent superposition of 0-quasi- 
particles with zero momenta, i.e. inflatons at rest. 

In this scenario, all the particles constituting the Uni- 
verse are created due to their interactions with the infla- 
tons. For a toy model describing the interaction between 
inflatons and other Bose particles Xi we may consider the 
term g 2 4> 2 x 2 01 9~ 2<7 4>X 2 m the Lagrangian. For Fermi 
particles ip we can include a Yukawa coupling htpffitp to 
the inflaton. If we assume for simplicity that the bare 
masses of the fields x an d ip are very small, the infla- 
tons will decay into ultra-relativistic particles. If the 
creation of these particles is sufficiently slow, the par- 
ticles will simultaneously interact with each other and 
come to a state of thermal equilibrium at the reheating 
temperature T r . The theory of this gradual reheating, 
developed long ago with the first models of inflation |Q , 
was based on the perturbative theory with respect to the 
coupling constants g 2 and h. Consider a quantum Bose 
field x with the eigenfunctions Xfe(t)e +lk x for comov- 
ing momenta k. The temporal part of the eigenfunction 



obeys an oscillator-like equation with a time-dependent 
frequency: tt 2 k (t) = j£ + g 2 <j)(t) 2 or g 2 <J<f>{t). The iter- 
ative solution to the equation of motion for Xk{t) gives 
the perturbation theory result for the particles' occupa- 
tion number: n& oc g <C 1. The resulting rate for the 
three-legs process 4> -> XX is given by T^ xx = 
By the same methods, one finds for the process <fi — > ipip a 
rate T<p^^ — ^ . Each of the created quanta carries 
away half of an inflaton's energy, This perturbative 
result can be interpreted as individual inflatons decaying 
independently of each other into pairs of x or tjj-'psxticYes. 

However, it was recognized that reheating can begin 
with a stage of boson production in a regime of paramet- 
ric resonance Indeed, in the equation for the eigen- 
mode Xk(t), inflatons act not as individual particles but 
as a coherently oscillating field 4>(t). The smallness of 
the coupling g 2 alone does not necessarily correspond to 
small occupation number nk- The oscillating effective 
frequency results in the parametric amplification 

of the modes: Xk ~ e^ kt . In this case, the energy of the 
inflatons is explosively transferred to the bose field. This 
process of bosonic preheating occurs very rapidly and far 
from thermal equilibrium. 

For fermions, the Pauli exclusion principle prohibits 
the occupation number from exceeding 1/2. For this rea- 
son, it has been assumed that fermions are created in 
the perturbative regime with the production rate 
and occupation number nk oc h 2 <C 1, as described in the 
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classical papers M . 

Let us, however, consider the Dirac equation for a 
massless quantum Fermi field t/j(t,x): 



[^V M - h(j>{t)\ t/> = 0. 



(1) 



Here, the inflatons producing fermions also act not as 
individual particles but as a coherently oscillating field 
4>(t). Therefore, we should check whether coherently act- 
ing inflatons produce the same number of fermions as an 
ensemble of individually decaying inflatons. It turns out 
that the actual production of fermions from the oscillat- 
ing inflaton field dramatically differs from the perturba- 
tive approximation. 

We begin our investigation with the model \\<t> A + 
h^><fn/j. A special feature of this theory is that the prob- 
lem of fcrmion production by the inflaton <f> in an ex- 
panding universe can be completely reduced to a similar 
problem in Minkowski space-time. Indeed, let us per- 
form a conformal transformation of the involved fields, 
ip = a4> and \P = a 3 / 2 ip, and use a conformal time vari- 
able, t = \/\{p 2 J ^|y. Then (without decay) the ampli- 
tude of the transformed field <p is constant, (p. The back- 
ground inflaton solution is given by an elliptic function 

ip(r) w (p cn (t, and the period of the background 

oscillations T ~ 7.416 in this model Q. 

The equation for the eigenfunctions of the quantum 
fluctuations in this theory can be reduced to a second- 
order equation for an auxiliary field X (r, x) , so that 

= [ij^Vfj, + hip] X. The eigenmodes of the auxiliary 
field have the form X k (T)e +l]sL ' x R r , where the R r are 
eigenvectors of the Dirac matrix 7 with eigenvalue +1. 
The temporal part of the eigenmode obeys an oscillator- 
like equation with a complex frequency which depends 
periodically on time 



Xh 



qf -iy/qf)X k = 



(2) 



The comoving momentum k enters the equation in the 



combination 



K 

\<p 2 



k 2 ; therefore, the natural units of 



momentum are y/Xip. The background oscillations enter 
in the form /(t) = cn (t, -^=\ having unit amplitude. 

The combination of the coupling parameters ^- = q ul- 
timately defines the solutions to Eq. @. 

Equation (Q) will be the master equation in our investi- 
gation. We choose the vacuum positive-frequency initial 
condition X/-(ro) = N k e~ lKT . The normalization factor is 

— 1 /2 

A4 = (2f2fe(f2fc + y/qfj) with the real part of the ef- 
fective frequency given by f2| = k 2 +qf(r) 2 . Notice that 
the problem of the production of bose particles \ m the 
theory j\(j) A + \g 2 4> 2 X 2 is reduced to an equation similar 
to (||) but without the imaginary part of the frequency 

and with the resonance parameter <? = @ • 

Using the standard ^-field operator expansion, we can 
express the comoving occupation number of particles in 



a given spin state through the solution of Eq. 
e.g., H for details) 



nk{T) = \-k Im 



Vqf 
2n k 



(see, 



(3) 



Certainly n k < 1/2 . The total comoving number den- 
sity of particles and antiparticles in both spin states is 
(2^)3 f d 3 k n k . Hence, the energy density of created 
fermions is = / d 3 k £l k n&. 

Eq. (||) can be solved numerically and the occupation 
number of particles can be found with (^). The results for 
^/c(t) with different parameters q are plotted in Fig. [|. 
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FIG. 1. The occupation number n k in A^ 4 -inflation as 
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a function of time r (in units of T) for q _ 
(lower), 1 (middle on right), and 100 (upper on right) and 
ft 2 = 0.18, 1.11, and 11.9, respectively. The period of the 
modulation (see Eq. (^)) is about 88, 20 and 22 (in units 
of T) accordingly. 

The occupation number exhibits high frequency (pe- 
riod < -j) oscillations which are modulated by a long 
period behavior. If we average the occupation num- 
ber over these high frequency oscillations, ftfc(r) = 

7> J^ T+T ^ d,Tnk(r), we can write the smoothed occupa- 
tion number of fermions in a factorized form 



rifc(r) = F k sin 2 v k r . 



(4) 



It is remarkable that, for all q, the occupation number 
of fermions evolves periodically with time || . The spec- 
trum (envelope) function F k and the period ^- depend 
on the parameter q. Typically ^- ^ 10 — 100, see Fig. ||. 
Therefore, the fermionic modes get fully exited around 
their time-average value ^ within just tens of inflaton 
oscillations! This is very different from the result of the 
perturbative approximation which gives the timing for 
excitation ~ T7_^^ which is about 10 14 inflaton oscilla- 
tions H . Notice that for q 1 , the number of fermions 
is almost constant between two successive zeros of the 
inflaton field. It jumps in a step-like manner at instances 
when (f>(t) crosses zero. This is a familiar picture for 
the creation of bosons explained in However, con- 

trary to the bosonic resonance which is unbounded, the 
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fermionic resonance reaches its peak value Ff. < 1/2 after 
a time and then falls back to zero, where it begins 
the cycle again. 

Fortunately, there is a semi-analytic shortcut for the 
calculation of Fk and v k . We found that the envelope 
function is given by the expression 



F, 



Here, X^(T) is the value of the first fundamental solu- 
tion of Eq. (||) (defined by the intial conditions xj^' (0) = 

1, (0) = 0) taken after one full oscillation; it is a func- 
tion of k which depends on the parameter q. The mod- 
ulation frequency v% is given by the relation cos^T = 
—ReX^ (T). Therefore, to find Fk and Vk, one need only 

calculate the complex value xjp (T) instead of perform- 
ing a full numerical integration of Eq. (Q) . 



1 



sin Vk 



T 2fif 



ImX$\T) 



(5) 




FIG. 3. The resonance excitation band of fermions in 
A0 4 -inflation for q = ^- = 1.0. The heavy line shows the 
envelope Fk of maximum occupation number calculated from 
Eq. (H). The light curve is the actual occupation number nk 
of each mode, calculated from Eq. (0), after 10 background 
oscillations. 
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FIG. 2. The envelope functions Fk showing the 

bands of fermion resonance excitation in A0 4 -inflation for 
q = K- = 10 _4 ,10 -2 , and 1.0 (the narrowest to broadest 
band, respectively). The band in the case q = 10 -2 already 
deviates considerably from the perturbative expectation. 

We calculated X^\t) numerically and constructed 
the envelope function Ff. plotted in Fig. |[ The spec- 
trum of nfc is very different from what is expected from 
the perturbative calculations, where it is narrowly peaked 
around k = y « 0.424 with the width r^.^. In Figs. §- 
^ we show, using (|), how the fermionic resonance bands 
are filled after 10 background oscillations. 

The function Vk gives us the time scale for fermion 
excitation. In Fig. @ we plot the period of modulation 
■jj- as a function of k. This function is peaked where 
Fk is peaked, i.e. the peaks of the resonance curve are 
the last to fill. For q > 1, the period for modes in the 
principle peak is about 20T. When q <C 1, modes in 
the principle peak grow more slowly, with a period of 
approximately . Interestingly, the period of modes out 
of the Ffc-peaks is significantly shorter. For this reason, 
the integrated number of fermions is quickly saturated to 
its time-average value. 




FIG. 4. As in Fig. 3, the resonance excitation band in 
A(/> 4 -inflation for q = h- — 100.0. The heavy line shows the 
envelope F^. The light curve is the actual occupation number 
of each mode rik after 10 background oscillations. 

The huge differences in the spectrum and rate of exci- 
tation of fermions produced by the coherently oscillating 
background field vs. those in the naive perturbative cal- 
culations are due to the interference of the fermion fluctu- 
ations. In this respect we can talk about the parametric 
excitation of fermions, where q plays the role of the res- 
onance parameter, and the frequency Vk can be viewed 
as the fermionic conterpart of the bosonic characteristic 
exponent fik 0. The amplitude of n% shows a distinct se- 
quence of resonance bands whose widths sharply narrow 
with increasing zone number. This is again reminiscent 



*Notice that, for the resonant bands in the bosonic case, the 
modulating factor in Eq. (^) is replaced by sinh 2 fi^t. This 
leads to the familiar estimation n*, ~ e 2Mfc ' for the number of 
bosons. 
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of the usual bosonic resonance. For q <C 1, the resonance 
bands are located at k = y (2/ + 1), where I = 0, 1, 2, 




FIG. 5. The log of the period of modulation 2- (in units 

of T) as a function of k 2 for q = ^ = 10 -4 , 1CT 2 , and 1.0 for 
the light, dotted, and heavy curves respectively. 

An important issue is the width of the parametric ex- 
citation of fermions. From Fig. || we notice that for any 
q the first, broadest zone has the width Ak ~ q 1 ^ 4 . 
This can be understood from the condition for non- 
adiabatic excitation of particles, Cl k > f2?, which gives 
us Ak < q^s/Xtp. This is similar to how the re- 
sult was elaborated for the broad bosonic resonance j|. 
There are also very narrow zones with higher momenta, 
which are also quickly saturated. This may be impor- 
tant for the production of superheavy fermions with a 
mass nrty ~ q^'^m^, where the effective inflaton mass 
~ 10 13 Gev. 

Let us estimate the energy that will be accumulated 
in created fermions. Averaging over short-period oscilla- 
tions, we have 

^ = ~2 J dkk 2 F k sin 2 v k T . (6) 

For the most interesting case of the broad resonance ex- 
citation, g > 1, we get ~ O.l/i 2 ^ 1 / 4 ^, where the 
inflaton energy is = jXf> 4 . In chaotic ^A<^ 4 -inflation, 
A ~ 10 -13 . If the resonance parameter q is large but 
the coupling parameter is small, h < 0.1, then only a 
small fraction of the inflaton energy will be transferred 
into fermions. However, if we push the parameters to 
h > 0.3 (q > 10 12 ), the inflaton energy can be transferred 
into fermions. This is quite an extreme case, because it 
requires SUSY cancellation of the radiative corrections 
~ h 4> . However, in the hybrid inflation scenario, for 
example, there is more freedom to chose the model pa- 
rameters which control the efficiency of preheating ||. 
There may be a range of parameters for which the infla- 
ton energy is fully converted into fermions even before 
the bosonic resonance develops. 

Let us now turn to the model ^m 2 (f> 2 4- h^p<j)tp, where 
m is the inflaton mass. In this inflationary model the 



inflaton evolution is simply described by harmonic os- 
cillations, <j) w $(t) sin mt, with decreasing amplitude, 
$(t) = - . However, the problem of fcrmion pro- 

duction in this model cannot be reduced to a problem 
in Minkowski space-time; the expansion of the universe 
is essential. Let us assume, for a moment, that there is 
no expansion and the amplitude $(t) = const. Then the 
master equation would be Eq. (^) where now the reso- 
nance parameter q = % and k 2 = (c.f. the bosonic 
resonance ||). The result for the occupation number of 
fermions in this case can also be found with Eq. (Q) and 
is plotted in Fig. 0. 
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FIG. 6. The resonance band in m 2 <j!> 2 -innation with 
q = h -n = 10 -4 , 10 -2 , and 1.0 for the narrowest to broadest 
band, respectively. The band in the case q = 10 already 
deviates considerably from the perturbative expectation. 

Without expansion of the universe, the qualitative re- 
sult is similar to that for the conformal theory. However, 
if we take the expansion into account, we can expect an 
adiabatic change in the structure of the fermion paramet- 
ric excitation. This is similar to the erasure of differences 
between resonant and non-resonant modes and the onset 
of stochastic resonance in the bosonic case ||. There- 
fore, we conjecture that, for the massive inflaton, the 
parametric excitation of fermions fills a Fermi-sphere of 
radius Afc ~ q 1 l 4 m without distinct resonance zones. 

If expansion redshifts the fermion modes fast enough, 
it can prevent them from being resonantly excited. We 
found that for qo < 1 or h < 10~ 6 we return to the 
perturbation theory result for the rate of fermion pair 
production ^r^(a 4n ^,^) ~ r^^^^n^, where the comov- 
ing number density of inflatons is = m$ 2 a 3 /2. This is 
to be interpreted as the creation of a pair of fermions, tp 
and ip, from the decay of a single inflaton at the instant 
of resonance between a given momentum k and inflaton 
mass m, k = ^ma(t). 

In this paper we studied the production of fermions 
by an oscillating inflaton field, and found that fermions 
are created in a regime of parametric excitation, which 
is very different from the perturbative regime. This can 
be important to many interesting cosmological applica- 
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tions including the channel of inflaton energy transfer, 
the estimation of the reheating temperature, the produc- 
tion of superheavy fermions, the production of gravitinos 
in the supergravity inflationary models |)|, etc. It also 
elucidates another feature of bosonic preheating, missed 
in the previous studies. Indeed, in the resonance bands 
bosons are exponentially unstable, rik ~ e 2Mfct . How- 
ever, outside of the resonance bands (Xk is imaginary and 
we will have bounded excitation of bosonic fluctuations 
rik ~ 0(1) given by Eq. (|J). How this will alter the 
preheating scenario and its lattice simulation, based on 
classical approximation of fluctuations (h]] , shall be seen. 
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